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The CDF detector simulation framework is integrated into an ACH — h application used to process events in the 
CDF experiment. The simulation framework is based on the GEANT3 package. It holds the detector element 
geometry descriptions, allows configuration of digitizers at run-time and manages the generated data. The 
design is based on generic programming which allows for easy extension of the simulation framework. The 
overall design, details of specific detector components and in particular the performance of the CDF simulation 
compared to collider data are described. 



1. Introduction 

The Collider Detector at Fermilab (CDF) Q op- 
erating at the Run II Tevatron Collider is a complex 
general purpose detector and has many subsystems. 
For the simulation of each subdetector, one needs ge- 
ometry, a digitizer to generate hits, a menu to change 
digitizer's configuration, and an object to keep a sim- 
ulated event. The simulation framework allows to 
plug these pieces into the existing software system and 
hides complicated details of their interaction from the 
user. To achieve these goals, the framework design 
is based upon a mixture of generic programming and 
object oriented principles. This mix, rather than a 
pure object oriented solution, allows the framework 
to be easily expandable while retaining as much time 
efficiency as possible. 



2. General overview 

The CDF detector simulation framework is inte- 
grated into an AC++ .2] application used to process 
events at the CDF experiment. The tracking of parti- 
cles through matter is performed by the GEANT3 
package. The CDF software uses the same geometry 
for event reconstruction and simulation. The CDF ge- 
ometry package provides the volume description and 
geometry tree creation. The CDF geometry was de- 
signed to allow seamless transition to GEANT4 if de- 
sired. It can be converted to a GEANT3 or GEANT4 
geometry. An output event of the simulation contains 
data in the same format as raw detector data plus 
Monte Carlo truth information. 

2.1. Functionality of cdfSim executable 

The main simulation executable, cdfSim, allows 
generation of physics events with different genera- 
tors such as Herwig v6.5, PYTHIA v6.2, Isajet v7.51, 
WGRAD, WBBGEN, GRAPPA (GRACE for pp), 
Vecbos, BGenerator, MinBiasGenerator, SingleParti- 
cle. Les Houches Accords - a universal interface be- 
tween matrix element generators and parton shower 



MC programs - is implemented in Herwig, PYTHIA 
and GRAPPA. 

In addition to decay routines internal to each MC 
generator, particle decays can be simulated by three 
decay packages - QQ v9.1, EvtGen and Tauola. 

The cdfSim executable may be used with differ- 
ent configurations of subdetectors, different geome- 
try levels, and physics processes depending on desired 
accuracy versus time efficiency. The following sub- 
detectors are simulated: 

• Silicon detectors (SVX, ISL) 

• Central Outer Tracker (COT) 

• Muon systems 

• Time-of-Flight system (ToF) 

• Calorimeters 

• Cherenkov Luminosity Counters (CLC) 

• Forward detectors (Miniplug, BSC, RPS) 

In addition, the CDF geometry contains a detailed 
description of passive material elements, in particular 
within the silicon detectors. 

2.2. Generators within AC++ framework 

The flow diagram in Figure shows how the var- 
ious generators are incorporated within the AC++ 
framework. The AC++ modules are shown in order 
of execution from top to bottom. A HEPEVT com- 
mon block is used for communication between differ- 
ent generators and decay packages. The Genlnput- 
Manager module creates an instance of the CdfRn 
class which interfaces to the CLHEP random num- 
ber generators. CDF unified the usage of the various 
random number generators used throughout the sim- 
ulation code including all generator packages. This 
ensures statistically independent production of large 
Monte Carlo samples on the CDF Production Farms. 
CdfRn manages the random engines of its clients and 
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Figure 1: Generator related modules within the ACH — h framework. 



allows restoring and saving of current engine config- 
urations. The GenlnputManager and RandomGen- 
Manager are followed by a sequence of generator mod- 
ules and decay modules. The combination of enabled 
modules from this sequence is checked at the AC++ 
framework level for user setup errors. Each module 
has an interface to the generator/decay package con- 
trolling the access to the underlying FORTRAN rou- 
tines. The generator or decay package modifies the 
HEPEVT common block. After the event has been 
processed by all enabled generators and decay pack- 
ages, the HEPEVT common block is converted into a 
persistent object and added to the event record. After 
this step, a sequence of generator level filters can be 
evoked to reject events based on user chosen selection 
criteria. Finally, the position of the primary interac- 
tion vertex is generated by the GenPrimVert module 
which allows to take beam line information from the 
CDF data base depending on a chosen run number. 
This information is also added to the event record. 

2.3. Simulation framework 

The diagram in Figure [21 shows the relationship be- 
tween different components of the simulation frame- 
work. The boxes represent ACH — h modules in order 
of execution from top to bottom. An instance of the 
CDF geometry is created by the GeometryManager 
module at the beginning of the job. The Simtnit- 
Manager module initializes GEANT3, passes the CDF 
geometry to GEANT3, and obtains the information 
which simulation configuration is used. The Simu- 



lationControl module instantiates the simulation el- 
ements requested by the SimlnitManager. The sim- 
ulation elements are processed by a simulation base 
which is implemented as an abstract factory owned 
by SimulationControl. Each simulation element con- 
sists of a geometry, a configuration menu, a digitizer, 
and an object to keep the Monte Carlo event data. 

During the event processing the SimulationControl 
module interacts with GEANT3 and passes to it the 
Monte Carlo particles from the generator. Inside an 
active detector volume, GEANT3 calls a user defined 
stepping routine and control is passed to the digitizer 
of the corresponding simulation element. This action 
is dispatched by the simulation base. The digitizer 
creates hits, depending on the configuration menu, 
and adds the information to the Monte Carlo event 
data. After tracing of particles with GEANT3 is fin- 
ished, the hit information from the MC event data ob- 
ject is converted into the persistent object and added 
to the event record. 

At the beginning of each simulation run, informa- 
tion from the CDF data base may be accessed to ob- 
tain noise, calibrations or alignment constants used 
for the simulation job. 



3. Simulation performance 

In the following subsections the current status of the 
simulation of the different subdetectors is described 
including a discussion of the simulation performance 
compared to collision data. 
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Figure 2: Generator related modules within the ACH — h framework. 



3.1. Silicon detector 

The Run II silicon detector consists of seven double 
sided layers and one single sided layer mounted on the 
beampipe covering a total radial area from 1.5-28 cm. 
The silicon vertex detector covers the full Tevatron lu- 
minous region and allows for standalone silicon track- 
ing up to a pseudo-rapidity ~ 2. The geometry 
of the silicon detectors has been implemented in great 
details in the simulation including passive materials 
and allows for detector misalignment. 

The main purpose of simulating the silicon detec- 
tor is to describe the charge deposition of traversing 
tracks on the silicon strips well enough to evaluate sil- 
icon tracking performance. The basic distributions to 
model are the number of strips included in a charge 
cluster and the cluster shape. 

Three charge deposition models (CDM) are imple- 
mented in the simulation framework: geometric, para- 
metric and physical. The geometric model is based 
purely on geometry. The total amount of charge de- 
posited on a sensor is calculated by GEANT3 based 
on an unrestricted Landau distribution. The amount 
of charge deposited on each strip is a fraction of the 
entire charge proportional to the path length of the 
particle in the vicinity of the strip. The result of this 
method is that for tracks at a certain incident an- 
gle, all strips in the cluster have the same amount 
of charge deposited. The cluster length and shape in 
silicon data is considerably more complicated due to 



various physics effects such as 5-rays. 

The parametric model is based on the physics ef- 
fects simulated by GEANT3. However, the GEANT3 
results have been parameterized into template distri- 
butions to increase the speed of the code. The im- 
portant effects that are included in this model are en- 
ergy deposition using a restricted Landau distribution 
with (5-ray production above a cutoff value, capacitive 
charge sharing between strips, and noise. This CDM 
has been further tuned to data. 

The physical model simulates the physics of charge 
deposition from first principals. The basic idea is to 
have the capability of monitoring, predicting and com- 
paring the electrical behavior of the real detector with 
the expectations from theory. The physics described 
by this model includes continuous energy loss using a 
restricted Landau distribution with (5-ray production, 
convolution of Landau fluctuations with a Gaussian, 
magnetic field effects, diffusion of electrons, holes, and 
noise. In addition to these effects, the charge inversion 
of a sensor can be simulated. 

In Figure |21 CDF silicon data are compared with 
simulation using the parametric CDM. The cluster 
length (left) and cluster profile (right) are compared 
for both the axial side (top) and stereo side (bottom) 
of Layer 1. The intrinsic silicon resolution is displayed 
in Figure 13 Unbiased axial residuals (Layer 2) are 
compared between data and MC. Silicon simulation 
and data are in good agreement. 
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Figure 3: Cluster length (left) and cluster profile (right) 
for the axial side (top) and stereo side (bottom) of 
Layer 1 of SVXII. Data are the black points, Monte 
Carlo are the shaded (red) points. 
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Figure 4: Unbiased axial residuals (Layer 2) . Date are 
the points with error bars; MC is the histogram. 



3.2. Central Outer Tracker 

CDF's Central Outer Tracker contains 30,200 sense 
wires arranged in 96 layers combined into four axial 
and four stereo superlayers. It also provides dE/dx in- 
formation for particle identification. The default drift 
model used in the COT simulation is based on the 
GARFIELD package [4J , a general drift chamber sim- 
ulation program. The default GARFIELD parameters 
are scaled to describe the data. 

An example of the performance of COT simulation 
is displayed in Figure |S] which shows a comparison 
of track residuals between data and MC. The track 
residual is calculated as the hit displacement minus 
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Figure 5: Track residuals (hit displacement minus track 
displacement) for data (points) and MC (histogram). 
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Figure 6: Signed track curvature (charge/pr) for primary 
tracks with impact parameter do < 3 mm in W — » [iv 
data (points) and MC events (histogram). 



track displacement after applying track quality cuts. 
The simulation (histogram) is in a good agreement 
with the data (points). 

Figure is a comparison of W — ► /iv events gen- 
erated with PYTHIA to W data obtained from an 
inclusive high pt muon data set with a muon trigger 
threshold of 18 GeV/c. Shown is the distribution of 
signed track curvature (charge/p-r) for primary tracks 
with impact parameter do < 3 mm. The data (points) 
are reasonably well described by the Monte Carlo (his- 
togram) . 

3.3. Time-of- Flight system 

One of the new devices in the CDF II detector is 
a Time-of-Flight system with a resolution of about 
100 ps. It employs 216 three-meter-long scintillator 
bars with photomultiplier tubes (PMT's) at each end 
located between the outer radius of the COT and the 
superconducting solenoid (see Figure 0) ■ The Time- 
of-Flight system is designed for the identification of 
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Figure 7: Schematic view of the CDF detector indicating 
the location of the scintillator bars and attached PMT's. 
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Figure 8: Comparison of time distribution between the 
data and simulation. Crosses with the error bars are 
data; the histogram is simulation. The upper (lower) 
plots are a comparison of the time distribution for the 
east (west) channels before (top) and after (bottom) 
time-walk correction. 



kaons with a 2 cr-separation between ir and K for p < 
1.6 GeV/c. 

The energy deposition of the incident particle as 
well as the time at entry and exit points are sim- 
ulated by GEANT3. This information is translated 
into TDC and ADC counts. Three response parame- 
terizations with different levels of detail in the param- 
eterization arc implemented. One of the models uses 
a simple analytic form to parametrize time and ADC 
count. Two other models use calibration parameters 
from ToF data. 

Figure |H1 shows a comparison of the time distri- 
bution between collision data and simulated events. 
Crosses with error bars denote data and the histogram 
represents the simulation. The upper (lower) plots 
are a comparison of the time distribution for the east 
(west) channels before (top) and after (bottom) time- 
walk correction. The comparison indicates reasonable 



Figure 9: Location of the CDF muon systems in azimuth 
<f> and pseudorapidity n. 

agreement between data and simulation allowing for 
further tuning of the simulation. 

3.4. Muon systems 

There are four muon sub-systems in CDF: Central 
Muon Detector (CMU), Central Muon Upgrade de- 
tector (CMP), Central Muon Extension (CMX), and 
Intermediate Muon Detector (IMU) . Figure shows 
the location of the muon systems in azimuth <f> and 
pseudorapidity 77. 

The central calorimeters act as a hadron absorber 
for the muon detection system. The CMU consisting 
of four layers of planar drift chambers is located di- 
rectly outside the central calorimeters. The CMU sys- 
tem covers \r)\ < 0.6 and can be reached by muons with 
Pt in excess of 1.4 GeV/c. To reduce the probability 
of misidentifying penetrating hadrons as muon candi- 
dates in the central detector region, four additional 
layers of drift chambers (CMP) are located behind 
0.6 m of steel outside the CMU system. To reach these 
two detectors, particles produced at the primary inter- 
action vertex must traverse material totaling 5.4 and 
8.4 pion interaction lengths, respectively. The CMX is 
located in the pseudorapidity interval 0.6 < |ry| < 1.0 
extending the polar acceptance of the muon system. 
This detector consists of free-standing conical sections 
at the four corners of the central detector with eight 15 
degree wedges. The IMU identifies muons over the ra- 
pidity range 1.0 < 1 77 1 < 1.5 by a 4-layer barrel of drift 
tubes (BMU) and scintillators (WSU, BSU, TSU). All 
muon systems have ADC/TDC read-out. The Run II 
upgrades to the muon system almost double the cen- 
tral muon coverage compared to Run I and extent it 
up to \r)\ ~ 1.5. 
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Figure 10: Top: Comparison of the number of CMUP 
muons for each stack in the North wall, for data (points 
with error bars) and the simulation (shaded histogram). 
Bottom: Ratio of the above distributions. 



The most challenging part in the muon simula- 
tion is the description of the complicated geometry 
of the muon systems, while tuning of the digitization 
to muon trigger data is straight forward. Special ef- 
forts were made to compare simulation to real data 
with the aim of checking the geometry description. 
The study was done for CMP, CMU, and CMX us- 
ing muons from W —* [iv and Z — > in data 
and simulation. The overall description of the data is 
good. As an example of the numerous studies, Fig- 
ure ^| shows a comparison of simulation and data for 
the number of muons per stack in the North wall. 



3.5. Cherenkov Luminosity Counters 

The Cherenkov Luminosity Counters p| are used to 
measure the Tevatron luminosity at CDF. The CLC 
acceptance for pp inelastic processes is estimated from 
the simulation and adds a major contribution to the 
luminosity uncertainty. An excellent CLC simulation 
performance is therefore critical. The generation and 
propagation of Cherenkov photons is fully simulated 
by GEANT3 and then corrected for photo-cathode ef- 
ficiencies. The hit count in the CLC is sensitive to 
the material traversed by particles between the in- 
teraction point and the CLC counters. Therefore, the 
detector geometry in front of and in the vicinity (back- 
scattering) of the CLC needs to be described with high 
accuracy. Monte Carlo and data are in good agree- 
ment as can be seen in Figure ITTl It shows a compar- 
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Figure 11: Amplitude distributions in the three layers of 
CLC counters. Dots are Monte Carlo, the red (west) and 
blue (east) histograms are data. 



ison of amplitude distributions in the three layers of 
the CLC counters. Dots are Monte Carlo (PYTHIA), 
the red (west) and blue (east) histograms are mini- 
mum bias data. The distributions are normalized to 
the single particle peak. The resulting uncertainty on 
the CDF luminosity calculation due to the CLC sim- 
ulation is less than 4%. 



3.6. Calorimeters 

The CDF calorimeters 0, 13 provide separate 
electromagnetic (EM) and hadronic (HAD) measure- 
ments as shown in Figure The central elec- 
tromagnetic (CEM) and hadronic (CHA) calorime- 
ters (| r\ | < 1.1) employ a projective tower geom- 
etry back to the nominal interaction point with a 
segmentation of Arj x A</> ~ 0.1 x 15°. The sam- 
pling medium is composed of scintillators layered 
with lead absorbers in the electromagnetic calorime- 
ter and steel in the CHA. The energy resolution 
for the CDF central calorimeter is (j(Et)/Et = 
[(13.5%/^/E^) 2 + (2%) 2 ] 1/2 for electromagnetic show- 
ers and a(E T )/E T = [(50%/V^) 2 + {3%) 2 } 1/2 for 
hadrons with Et measured in GeV. A layer of pro- 
portional chambers (CES), with wire and strip read- 
out, is located six radiation lengths deep in the CEM 
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Figure 12: Elevation view of one half of the CDF detector 
displaying the various calorimeter compartments. 



\n\ Range 


A0 


An 


0. - 1.1 (1.2 h) 


15° 


~ 0.1 


1.1 (1.2 h) - 1.8 


7.5° 


~ 0.1 


1.8 - 2.1 


7.5° 


~ 0.16 


2.1 - 3.64 


15° 


0.2 - 0.6 



Table I CDF Run II calorimeter segmentation. 



calorimeter, approximately near shower maximum for 
electromagnetic showers. The CES provides a mea- 
surement of electromagnetic shower profiles in both 
the 0- and z-directions. Proportional chambers lo- 
cated between the solenoid and the CEM comprise 
the central preradiator detector (CPR), which samples 
the early development of electromagnetic showers in 
the material of the solenoid coil, providing informa- 
tion in r-0 only. In addition, the forward calorime- 
ters have been replaced in Run II by a new scintilla- 
tor tile based plug calorimeter which covers the range 
\rj\ < 3.6 and gives good electron identification up to 
1 77 1 ~ 2. The plug electromagnetic calorimeter also 
has fine grained shower profile detectors at electron 
shower maximum, and preshower pulse height detec- 
tors at approximately IX D depth. The calorimeter 
cell segmentation is summarized in Table [I] A com- 
parison of the central and plug calorimeters is given 
in Table HD 

The main objective for the calorimeter simulation 
is speed. The simulation of the CDF calorimeters 
is based on the GFLASH 9] parameterization pack- 
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Table II Comparison of CDF Run II central and plug 
calorimeters. 



age interfaced with the GEANT3 simulation. The 
simulation of EM and HAD showers in GFLASH 
is initiated when particles undergo inelastic interac- 
tions inside the calorimeter volume. GFLASH uses a 
mixture-level GEANT3 geometry description for the 
CDF calorimeters. It handles the spatial distribution 
of the deposited energy and the energy loss of particles 
inside the geometry volumes. 

The GFLASH parameters are tuned. The basic idea 
of the tuning process is first to bring the response of 
the electromagnetic calorimeter into agreement with 
data and then tune the hadron response by adjusting 
GFLASH parameters to match pion tcstbeam data 
and pp collision data. The tuning of GFLASH is split 
according to the two sets of parameters that control 
the fraction of visible energy produced in the active 
medium, and the energy and depth dependent spatial 
distribution of the various components of the shower 
model. 

To tune EM showers, 8-227 GeV electrons from test 
beam data were used. The energy scale and resolution 
for electrons was verified using E/p. The hadronic 
shower shape was tuned to test beam data and min- 
imum bias data. High pt tuning was done using 
high pt pions test beam data in the energy range 7- 
227 GeV (mainly 57 GeV) for CHA and 8-227 GeV 
for PHA. The low pr tuning process uses isolated low 
Pt tracks from minimum bias data. The Minimum- 
Ionizing-Particle (MIP) peak is measured in the EM 
compartments with pions. The hadronic energy scale 
in the tuning procedure is set by the response of 
57 GeV pions in the central and plug calorimeters. 
Since there are no test beam data available for the 
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Figure 13: Comparison of E/p distributions between 
simulation and testbeam as well as collision data for 
different particle momenta. 
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Figure 14: Comparison of energy depositions in the CDF 
hadronic (left) and electromagnetic (right) calorimeter 
for 4 GeV/c (top) and 40 GeV/c (bottom) muons. 



wall calorimeter, the energy scale of WHA is set to 
CHA. Once the MIP distribution and the hadron en- 
ergy scale are set, the logarithmic energy dependence 
of the calorimeter response can be tuned. Finally, 
for the tuning of the lateral shower profile, tracks ob- 
tained from minimum bias events measured in the cen- 
tral part of the CDF detector are used. The available 
energy range was 0.5-2.5 GeV. 

Figure [21 Figure ^3 shows a comparison of the 
E/p distributions between simulation and testbeam 
as well as collision data at low p for different parti- 
cle momenta. The calorimeter response in the simu- 
lation yields excellent agreement with data. As an- 
other example of calorimeter simulation performance, 
compares the energy deposition of muons in the CDF 
calorimeter. Data are shown as dashed lines and MC 
as solid. The two top plots show the energy distribu- 
tions for 4 GeV/c muons from J/tp — > fJ, (J>~ in the 
hadronic (left) and electromagnetic (right) calorime- 
ters. The two bottom plots show these distributions 
for 40 GeV/c muons from W — > fiv decays. Again, 
excellent agreement is seen between data and Monte 
Carlo. 

The GFLASH based calorimeter simulation is about 
one hundred times faster then the GEANT3 shower 
simulation for a typical top quark event. 



4. Conclusion 

We described the overall design of the CDF de- 
tector simulation framework which is based on the 
GEANT3 package and integrated into an AC++ ap- 



plication used to process events in the CDF experi- 
ment. The CDF simulation framework is shown to be 
flexible, easily extensible, and efficient. It is hiding 
complex infrastructure details from a user and suit- 
able for large scale Monte Carlo production. Subde- 
tectors are successfully implemented within the CDF 
simulation framework. We discussed details of the 
simulation of specific detector components and in par- 
ticular the performance of the CDF simulation which 
shows good agreement with pp collider data. A fur- 
ther upgrade to GEANT4 is possible through the low 
level API without affecting the client's code. 
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